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We study a class of supersymmetric models with non-universal gaugino masses that could 
arise from _F-terms in a general combination of the singlet and adjoint representations of 
SU(b). We explore models that satisfy present Large Hadron Collider and other bounds, 
showing how the allowed parameter space is divided into distinct "continents" . Regions of 
parameter space that ameliorate the supersymmetric little hierarchy problem with a small 
fi parameter include the usual focus point scenario, but also natural areas with much lighter 
squarks and sleptons. These models are continuously connected in parameter space to regions 
in which stau co-annihilation or Higgs exchange is mostly responsible for dark matter annihi- 
lation, and to models in which the thermal relic abundance is achieved by slepton-mediated 
annihilation, reviving the bulk region that is severely restricted in mSUGRA models. In 
hybrid or confluence regions, several mechanisms combine to give the requisite dark matter 
annihilation rate. In each case we study the prospects for direct detection of dark matter. 
We also comment briefly on the impact of recent hints for near 125 GeV from the LHC. 
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I. INTRODUCTION 

Super symmetry [l| as an extension of the Standard Model (SM) provides a way of explaining 
why the weak scale is so small compared to the Planck scale or other very high energy scales in 
fundamental physics. The Large Hadron Collider (LHC) is presently engaged in searches with 
increasing sensitivity to supersymmetry 0, 0] , excluding significant areas of parameter space cor- 
responding to lower gluino and up-squark and down-squark masses. Searches by underground 
detectors for a stable neutralino lightest supersymmetric particle (LSP) as the dark matter are 
also probing interesting parts of parameter space 0, @|. In addition, the non-observation of a 
Higgs boson at the CERN LEP e + e~ collider 0] places a very non-trivial constraint on the su- 
persymmetric models. These searches increase the lower bounds on superpartner masses, leading 
to an apparent fine-tuning problem known as the supersymmetric little hierarchy problem. The 
essence of this problem is the tension between the Z boson squared mass and the much larger soft 
supersymmetry-breaking squared masses. 

In the Minimal Supersymmetric Standard Model (MSSM), this problem can be understood by 
considering the relationship between input parameters and the electroweak scale, which may be 
written as 

m| = -2(|^| 2 + m 2 H J + ... (1.1) 

where the ellipses denote loop correction effects (which can be made small by an appropriate choice 
of renormalization scale) and a tree- level contribution suppressed by 1/ tan 2 /3. Here tan/3 is equal 
to the ratio of Higgs expectation values {H u } / {H^} , /i is the superpotential Higgs mass parameter, 
and mjj is the soft supersymmetry breaking squared mass of the Higgs boson that couples to 
the top quark. The cancellation needed between |/x| 2 and may therefore be regarded as an 
indication ^ of the fine-tuning required to obtain the observed weak scale. 

In terms of the input running gaugino masses Mi, Mi, M3 at the apparent unification scale 
M\j = 2 x 10 16 GeV, one finds from running two-loop renormalization group equations [7] that 
at the TeV scale is approximately: 



m 



\ u = 1.82MI-0.21MI+0.16M3M2 + 0.023M1M3 + 0.006M1M2 - 0.006M! 2 

-0.32A)M 3 - 0.07A M 2 - 0.022m 2 ,. (1.2) 

Here we have chosen tan (3 = 10 and chosen flavor-blind scalar squared masses m 2 , and a universal 
scalar cubic coupling parameter Aq for illustration. In the subspace of models with unified gaugino 
masses M 1 = M 2 = M 3 = m 1/2 at M v (often referred to as either "mSUGRA" or "CMSSM" in 
the literature), lower bounds on the gluino mass imply a significant fine tuning in order to reconcile 
equations (II. 1|) and (|1.2p . This will worsen if the LHC sets new limits, a plausible prospect in the 
very near future. 

This motivates models with non-universal gaugino masses, in particular those in which the 
gluino mass parameter M3 is relatively small compared to the wino mass parameter M 2 . As can be 



' Because there is no objective measure on parameter space, we do not attempt any detailed quantification of 
fine-tuning, but do view it as a qualitatively valid motivation and concern. 
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seen from eq. (|1.2|) . a small ratio of M3/M2 can lower —m? H and therefore decrease the amount of 
cancellation needed with |^| 2 [8]. In this paper, we will explore some features, including dark matter 
properties, of a class of non-universal gaugino mass models, with particular attention to models 
that can naturally accommodate small and are therefore more attractive from the perspective 
of the super symmetric little hierarchy problem. 

Specifically, we will use a class of models that have the gauginos getting contributions to their 
masses from an F-term that transforms as a 24 representation of the global SU (5) group that 
contains the SM gauge group [{J]. This yields a contribution to the gaugino masses in the ratio 
Mi : M2 ■ M% = 1 : 3 : — 2. The same pattern arises if the F-terms transform as a 54 
representation of the global 50(10) group that contains SU(5). In either case, the group may or 
may not be promoted to a true grand unification gauge group. Gaugino masses in this pattern 
may be added to a universal contribution M\ : M% : M3 = 1:1:1. 

In the CMSSM, there are four main mechanisms that allow for dark matter annihilation that 
is efficient enough to prevent the early universe from becoming matter dominated too soon. First, 
the "bulk" region at small mo and small mi/2 allows efficient dark matter annihilation by t-channel 
slepton exchange. This region is under considerable pressure (if not eliminated entirely) by the LEP 
and LHC bounds, but with gaugino mass non-universality we will see that it is easily resurrected. 
Second, there is a co-annihilation region at small mo, in which the lightest stau co-exists and co- 
annihilates in the early universe with an LSP that is not much lighter [ll|. Third, a "focus point" 



or "small fj," region 12| at very large mo, has \i small, so that mostly bino-like LSPs contain a 
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subdominant but significant Higgsino admixture, allowing them to annihilated efficiently 
and through weak bosons. If the LSP mass exceeds mt, then the final state of dark matter pair 
annihilations can be tt in this case. Fourth, at large tan fi the ^4°-funnel region allows the LSPs to 



annihilate through s-channel pseudo-scalar Higgs exchange 13j. In the more general MSSM, there 
is also a possibility of LSPs that co-annihilate efficiently because they are wino-like [14j, but this 
cannot be realized in the unified gaugino mass case. In the following, we will take the dark matter 
density to be in the conservative range 

0.09 < n BM h 2 < 0.13, (1.3) 

where h ~ 0.71 is the Hubble parameter today in units of km/(sec Mpc). It is very important 
that it is optional to require the predicted thermal relic abundance to lie in this range, because 



there are a variety of easy ways to evade this bound (see for example [15|). If the predicted f^DM^ 2 
for a particular model comes out lower than 0.09, then axions or some other species could make 
up the difference, and so this region of parameter space should not be viewed as disallowed. If 
SIdm^ 2 > 0.13, then one can reduce the predicted density by having the apparent LSP decay to 
a lighter singlet particle, or by invoking imparity violation so that there is no supersymmetric 
dark matter at all. Still, it is interesting to consider the rang e eq . (jl.3p as providing a minimal 
accommodation of astrophysical and cosmological observations 16]. 

Modifying the universal boundary conditions of CMSSM can lead to interesting new conditions 
for dark matter; see for some recent examples [l3]-[2^] and references therein. In [3(J, it was 
suggested that with a mixture of the two gaugino mass patterns mentioned above, chosen to 
make [i smaller, the LSP remains mostly bino-like but the observed thermal relic abundance of 
dark matter may be naturally obtained by efficient annihilations N\N\ — > tt through t-channel 
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exchange of top squarks. In these "compressed super symmetry" models, the ratio of the masses of 
the heaviest and lightest superpartners is reduced and \fx\ is also much smaller than encountered in 
the CMSSM. However, only models continuously connected in parameter space to the CMSSM case 
were considered in 30f | . In this paper, we will study the full parameter space obtained by varying 
over the entire allowed range spanned by linear combinations of the two gaugino mass contribution 
patterns M\ : M% : M3 = 1:1:1 and 1 : 3 : — 2. The viable parameter space is typically divided 
into three separate "continents", with distinctive features. In the following, we will characterize 
these regions in terms of the size of the fi parameter, and in terms of the mechanisms that can 
allow the dark matter abundance to be in approximate agreement with cosmological observations. 
In the larger parameter space spanned by non-zero F terms in the 24 of SU(5), the dark-matter 
allowed regions mentioned above merge and are deformed in interesting ways. One of the more 
interesting scenarios is a parameter-space region with small that is continuously connect to 
the CMSSM focus-point region, but occurs at much lower values of mo when gaugino-mass non- 
universality is significant. There are similar regions that are not continuously connected to the 
CMSSM focus point region. There are also stau co-annihilation and bulk regions, some of which are 
not continuously connected to the corresponding CMSSM regions, and which have very different 
prospects for dark matter direct detection. We will also discuss an intriguing viable region formed 
by a confluence between the A°-funnel, small-//, and stau co-annihilation regions, which can occupy 
a large chunk of parameter space, and which we refer to as the confluence island. However, this 
region is severely impacted by indirect constraints, as we will see. 



II. PARAMETERIZATIONS AND CONSTRAINTS 

In this paper, we will parameterize a general combination of the universal and non- universal 
gaugino mass patterns mentioned in the Introduction in terms of an overall gaugino mass scale 
mi/2 and an angle 024, which we define as: 

Mi = m 1 / 2 (cos #24 + sin 2 a) , (2.1) 
M 2 = m 1/2 (cos 024 + 3 sin 2 4), (2.2) 
M3 = mi/2 (cos 024 — 2 sin 6*24) . (2.3) 

Note that 024 = corresponds to the usual unified gaugino mass scenario, while 024 = ±7r/2 
correspond to a pure 24 of SU(5) [or 54 of 50(10)] F term contribution. Taking 024 024 + tt 
flips the signs of all three gaugino masses simultaneously, which is physically the same as the 
original model with the signs of the scalar cubic couplings and (i term flipped. 

Because of the high dimensionality of parameter space, it is not possible to do a complete study, 
and some choices must always be made in order to illustrate qualitative features. In this paper, 
we choose to include only a universal scalar squared mass mfj, despite the fact that non- universal 
gaugino masses might be expected to be accompanied by non-universality of some kind in the 
scalar sector. We assume this scalar mass non-universality to be small, or that it would lead 
to qualitative features similar to the results we obtain, but it would certainly be interesting to 
consider alternatives. In this regard, note that the scalar soft squared masses will arise in the form 
F*F/Mp^ nck , which transform in the product of the representation Rp of F and its conjugate, 
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Rp. The direct product representation Rp X R* F will always include the SU(5) [or SO(W)] singlet 
representation, so that it is sensible to take a common tbq, but it also can include other non-trivial 
GUT representations. Therefore, it would be just as reasonable to explore other patterns, but we 
will not do so here. The parameters mo and 024,/tt will be varied independently. 



In our explorations of parameter space, we used SOFTSUSY 3.1.2 [3l( to generate supersymmet- 
ric spectra,* using rru = 173.3 GeV, M b (M B ) MS = 4.25 GeV, and a s {M z ) = 0.118. The program 
micrOMEGAs 2.2 [331 ] was used to evaluate dark matter thermal relic abundance, the LSP-nucleon 
cross-section for dark matter detection, and implement constraints from other observables as fol- 



lows. For B — y tv mediated by charged Higgs bosons, we follow [341 ] by taking a constraint 

m H + > 13.5 tan /3. (2.4) 

We also consider a constraint 

BR(B S -> /x + aO < 11 x 10~ 8 (2.5) 

from [35J. There is a well-known deviation of the anomalous magnetic moment of the muon 



measured by |36[ from the Standard Model predicted values. Since we are unwilling to interpret 
this measurement as ruling out the Standard Model, we follow the "super-conservative" attitude 
of 37[ by only requiring agreement within 5<r: 

- 2.5 x 1(T 9 < Aa^ < 6.5 x 1(T 9 (2.6) 

Here Aa^ is the super symmetric contribution to (g^ — 2)/2. In the models we consider, this only 
impacts the very high tan /3 case. For b sj, we implement a constraint 

2.0 x 10~ 4 < BR(6 -> s7) < 5.0 x 10" 4 , (2.7) 



which is also somewhat more conservative than commonly imposed [38]. This is because constraints 
on the supersymmetric contribution to this observable only apply if one accepts the additional 
assumption of minimal flavor violation in the soft super symmetry-breaking sector. Because there 
is no reason not to expect additional sources of flavor violation, this constraint should not be 
considered mandatory. One could even adopt the point of view that no strict BR(6 — > sj) limit 
applies at all, although relying on non-minimal flavor violation to counteract the minimal flavor 
violating contribution might be seen as requiring some fortuitous tuning. In the models encountered 
below, h° has couplings very close to those of a Standard Model Higgs boson, but we use a slightly 
lower bound than the Standard Model 114.4 GeV LEP bound because of theoretical uncertainties 



In some cases, especially at high tan/3 and large mo, we found significant numerical differences between the model 
calculator we used, SOFTSUSY, and an alternative, SuSpect [33], evidently due to high sensitivity of the spectrum 
to the particular treatment of radiative corrections and translation of input parameters into physical masses. 
However, we expect the qualitative features to remain robust. 
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in the mass prediction in supersymmetry; we take 

m h o > 113 GeV. (2.8) 

All of the constraints mentioned above are implemented as contours in plots, rather than by re- 
moving the models from consideration, so that the impact of these indirect observables on the 
parameter space can be judged by the reader. Finally, we note that the CMS and ATLAS collab- 
orations have announced significant limits 0, S] on gluino and squark production, although only 
in the context of unified gaugino mass models and other simplified models with large hierarchies 
between the gluino and LSP masses. Most of the models we will study evade the LHC limits by 
taking sufficiently heavy gluinos, with M3 held fixed. In cases where Mi is held fixed instead, we 
will also plot the contours for rriq = 800 and 1000 GeV for the average squark mass on graphs 
below, as a very rough indication of where the exclusion reach of present LHC data might lie. 
However, because the superpartner mass spectrum is compressed in many of the parameter space 



regions we study, and the reach of the LHC can be significantly reduced in such cases [39(], a true 
exclusion would require a dedicated study that is beyond the scope of this paper. 



III. EXPLORATIONS WITH FIXED M 3 

In this section, we study several slices of parameter space with M3 held fixed at 600 GeV, so 
that the gluino mass is beyond the present reach of the LHC, in the range 1350 < M g < 1525 



depending on the masses of squarks (which contribute to the gluino mass in loop corrections [401]). 
We will consider both moderate and large tan/3, and varying mo and #24- From eqs. (I2.ip -( [2T3]) . 
one has: 

i\j */r ( ! + t an6l 24\ „, j, r A + 3tan6> 24 \ ,„ ^ 

Ml = M M 1 — a~ J ' M 2 = M 3 — — — . (3.1) 

\ 1 — 2 tan #24 / \ 1 — 2 tan #24 / 

Since we are holding M3 fixed in this section, \M\\ becomes very small when 024 /tt approaches 
— 1/4 and 3/4. Therefore, in the graphs below we take the range for 024,/tt from —1/4 to 3/4; 
this range (rather than, say, to 1, or —1/2 to 1/2) avoids splitting up the continents of viable 
parameter space in unnatural ways. At 624/^ ~ —0.102, one has M2 approaching 0, so models 



near this vertical line in figures below are excluded by LEP2 limits on charginos 38], providing 
an ocean between the two continents. At 624/11 ~ 0.148, one has \M\ \ and IM2I diverging without 
bound, so models near this vertical line are likewise excluded, providing another gap between the 
two continents of viable parameter space on either side. The boundaries of the continents on either 
side of this line are actually determined by the fact that small \M$/Mi\ leads to small |/i| 2 , as can 
be seen from eqs. (jl.ip and (|l,2p . 

This is illustrated in Figure HJ which maps the viable regions of parameter space in the (#24, fno) 
plane, for fixed M3 = — Aq = 600 GeV, tan/3 = 10, and \i > 0. The choices of Aq and tan/3 are 
motivated by the need to obtain a lightest Higgs mass above the LEP bound, and indeed all of 
the indirect bounds eq. (j2.4j) - (j2.8|) are satisfied throughout the shaded regions shown. Regions 
left blank do not have viable electroweak symmetry breaking, do not have a neutralino as the 



LSP, or have charged superpartners (staus or charginos) that are below the LEP2 bounds 38] 
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FIG. 1: Maps of the /i parameter (upper plot) and the predicted thermal relic abundance of dark matter 
Odm^ 2 (lower plot) in the #24, mo parameter space with fixed A/3 = — Aq = 600 GeV, tan j3 = 10, and fi > 0. 
In the upper plot, the dark (black) regions at the top and closest to the center represents /1 < 300 GeV, and 
successive regions lower and away from the center correspond to 300 GeV < fj, < 400 GeV (brown) and 400 
GeV < fi < 500 GeV (red) and so on up to the last region (blue) which corresponds to 800 GeV < \x < 900 
GeV. The regions left blank do not have viable electroweak symmetry breaking, do not have a neutralino as 
the LSP, or have superpartners that are too light, as described in the text. In the lower plot, the thin dark 
regions (black) correspond to the observed range 0.09 < Udm^ 2 < 0.13. The large interior regions (yellow) 
correspond to f^DM^ 2 > 0.13, while the darker shaded exterior region (green) has fiDM^ 2 < 0.09. 

The lower boundaries of the allowed continents in Figure [TJ and similar figures below, are set by 
the requirement that the LSP is a neutralino and not a stau, as required for a super symmetric 
explanation of the dark matter. The upper boundaries come from the requirement of proper 
electroweak symmetry breaking with positive |/u| 2 in the scalar potential. On the outside parts of 
the plots, -0.25 < 6» 2 4/V < -0.16 and 0.66 < 6» 24 < 0.75, the LSP mass will be below m z /2. 

The top plot in Figured] is a map of the /i parameter, with lower values arguably corresponding 
to less fine-tuning, as described in the Introduction. Note that the MSSM corresponds to the 
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vertical line #24 = 0. On that line, the smallest values of fi occur at the largest allowed values of 
the scalar masses, near mo = 4300 GeV; this is the focus point solution. For smaller values of mo, 
the solution for electroweak symmetry breaking gives larger values of fi, but it never exceeds 900 
GeV throughout the whole plot. The focus point region is continuously connected to a small-/x 
region near #24/ V = 0.65 with a wide range of mo values. This region connects to a region where 
Qdm^ 2 < 0.09 which extends down to very low values of mo < 100 GeV. These models survive 
the requirement of a neutral (non-stau) LSP because here [i < Mx,M%, so that the LSP is a 
higgsino-like state that is much lighter than the sleptons. 

The left continent in Figure [1] corresponds to M2/M3 < 0, and contains small- ^ regions only 
at large values of mo, comparable to the CMSSM focus point. Recall that the far left of the plot 
corresponds to very light bino-like neutralino LSPs, and that the region is separated from the 
central, CMSSM-like continent by a region where M2 is too small, leading to charginos below the 
LEP bound. The boundaries of that region are very nearly vertical, because the chargino mass 
depends only very weakly on mo through loop effects. 

The right continent in Figure Q] has both M1/M3 and M2/M3 negative, and is separated from 
the central CMSSM-like continent by a region where IM2/M3I is very large, as discussed above. In 
fact, the region from 0.08 < 624/ 'vr < 0.33 is excluded because — rn 2 H ^ is negative at the TeV scale, 
precluding electroweak symmetry breaking. The top edge of this continent is again a small- /j, focus 
point-like region, and the left edge of the continent near 624/ n = 0.35 also has small over a large 
range of mo- In fact, one can see from the plot that the range of parameters over which \x is small 
is considerably broader than the corresponding region on the central CMSSM-like continent, so 
that arguably this region is less fine-tuned in producing correct electroweak symmetry breaking. 

The lower plot of Figure [U shows the same regions, but now indicating whether the predicted 
thermal relic abundance of dark matter is larger, smaller, or within the range of eq. (|1.3p . The 
region where eq. (jl.3j) is satisfied is quite thin; this does not reflect any fine-tuning, but rather 
the fact that the dark matter relic abundance is now quite precisely known from experiment. In 
the central continent, the interior predicts too much dark matter, due to inefficient annihilation, 
while on the boundary there is too little dark matter. For the CMSSM at #24 = 0, the two allowed 
possibilities are the focus point region near mo = 4325 GeV and the stau co-annihilation region 
near mo = 140 GeV. The top and bottom edges of the central continent share these features. The 
vertical left edge of this continent has efficient wino co-annihilations; because of small M2 the LSP 
has a significant wino content. The right edge achieves the correct dark matter abundance due to 
having a significant higgsino content of the LSP because [i is small, similar to the focus point case. 
Note that these regions are actually continuously connected, with hybrid features at the corners. 

In the right continent of Figure [TJ the top and left edges again have small ^. The bottom 
region is stau-coannihilation up to about #24 = 0.58, where it continuously merges into a fatter 
"bulk region", characterized not by stau coannihilations but by slepton-mediated annihilations to 
leptons. Note that the stau-coannihilation region has [i between about 300 GeV and 700 GeV, while 
the bulk region has \i between 700 and 800 GeV. Also visible as two small peaks are /i°-resonance 



411 ] and Z° resonance regions near 624/n = —0.15 and —0.17 respectively. The bulk region points 



at small mo on the right have light slepton masses not far above their LEP2 bounds. For the 
left continent of Figure [H the h° resonance region cuts vertically through the entire continent, 
connecting to small-// and bulk regions at the top and bottom of the continent, respectively. Also 
shown on Figure [1] are eight selected example models labeled A, B, C, ... H, for future reference. 
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FIG. 2: The ratio of the spin-independent LSP-nucleon cross section to the current XENON100 limit for 
model points in Figure Q] (with fixed M3 = — Aq = 600 GeV, tan/3 = 10, and fi > and varying 824, mo) 
that have rioM^ 2 = 0.11. Different symbols are used for the model points according to which dark matter 
annihilation channels are most important in the early universe. 



For the models that do lie in the range of eq. (jl.3p for Qdm^ 2 ) it is interesting to consider the 
searches for direct dark matter detection, with the most sensitive at this writing being the results 
from XENON100 in ref. Q, which superseded earlier results from CDMS 0]. In Figure [fl we 
show the ratio of the spin-independent LSP-nucleon cross section to the current XENON100 limit 
for that LSP mass, for model points in Figure Q] with Odm h 2 = 0.11. The models are denoted by 
different symbols depending on what is the most important mechanism for dark matter annihilation 
in the early universe, and plotted both as a function of fj, and of 624- However, it is important to 
realize that the models with ct^z/ctxenonioo > 1 are not ruled out, for several reasons. First, the 
would-be LSP may have decayed, either to a lighter singlet or by i?-parity violation, or may have 
been diluted by late inflation. Second, there are significant uncertainties in the local density of 
dark matter, and especially in the nuclear matrix elements that go into the computation of asi- 
These could easily reduce the true &si by more than a factor of 2. Therefore, none of the models 
considered can be taken to be ruled out, but we see that within a default interpretation, most of the 
small-// models in the CMSSM continent, and some of those within the right continent can be said 
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FIG. 3: Comparison of superpartner and Higgs mass spectra for four selected model points A,B,C,D (top 
left, top right, bottom left, and bottom right, respectively) in the small- fi dark matter regions in Figure [TJ 



to be challenged by the XENON100 results. In contrast, the stau co-annihilation, bulk region, and 
wino co-annihilation models are not at all challenged by the present direct dark matter searches 
in this model set. In the case of stau co-annihilation models with #24/V > 0.35, asi does not fall 
monotonically with fi, but instead falls sharply and by orders of magnitude where it is subject to 
accidental cancellations between the light and heavy Higgs contributions and other contributions. 
These cancellations are features of a leading order calculation, and including loop corrections and 
real emission diagrams should be expected to remove these accidental cancellations, but osi can 
still be so small that direct detection will be quite problematic. The eight example model points 
from Figure [Dare labeled individually. Note that of the small-^ models, those with 624/^ < 0.45 
most easily evade the dark matter direct detection searches. 

As examples, we show in Figure [3] the superpartner and Higgs mass spectra for the four selected 
models with small fi and Odm^ 2 within the WMAP range, as labeled in Figures [1] and [2j The 
parameters for these models are M3 = — Aq = 600 GeV, tan/3 = 10, and: 



A: (^/t, mo) = (0.0, 4325 GeV) focus point, central continent (3-2) 
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B : (#24/71") mo) = (0.0671, 400 GeV) small //, central continent (3.3) 
C : (024/tt , mo) = (0.346, 400 GeV) small /i, right continent (3.4) 
D: (#24/tt, m ) = (0.45, 3440 GeV) focus point, right continent (3.5) 



Model A is a CMSSM focus point model, with the well-known features of squarks and sleptons 
that are far too heavy to be seen at LHC, higgsino-like neutralinos and charginos that are lighter 
than their wino-like counterparts, and a Higgs boson that easily evades LEP bounds. Model D is 
qualitatively similar, although with a slightly lighter spectrum. This implies that in the case of 
model A, the final state of dark matter annihilation is mostly tt, while it is mostly WW and ZZ 
in model D. In both of these models, discovery at the LHC will eventually come from gluino pair 
production to multi-jet final states. 

In contrast, models B and C both have squarks that will be accessible to the LHC. Both models 
have significant hierarchies between the right-handed and left-handed sleptons, and between the 
right-handed and left-handed squarks. Eventual measurements of these masses would help to 
confirm the role played by large M2/M3 in these models. Both models also have small mass 
differences between the higgsino-like neutralino and chargino states N2, N$, C\ and the LSP. In the 
case of model B, co-annihilations of all of these states are important in the early universe, while in 
model B the annihilation is mostly LSP pairs to it and ZZ. We also note that the spin-independent 
nucleon-LSP cross-section is a factor of 4 smaller for model C than for model B, and the former is 
below the limit from XENON100 at that mass while the latter is above the limit, when computed 
using the default micrOMEGAs parameters. However, another key difference between these models 
is that the lightest Higgs mass rrih is in tension (near or possibly below) the LEP limit for model 
C, while it is about 119 GeV for model B. 

We also show in Figured] the superpartner and Higgs mass spectra for the four selected models 
E,F,G,H with small mo, leading to stau co-annihilation and slepton-mediated dark matter regions 
within the Odm h 2 range from WMAP. The parameters for these models are M3 = — Aq = 600 
GeV, tan (5 = 10, and: 



E 
F 
G 
H 



(#24/tt, m o) = (0.0, 140 GeV) stau co-annihilation, central continent (3.6) 

(#24/71", mo) = (0.045, 215 GeV) stau co-annihilation, central continent (3.7) 

(#24/71", mo) = (0.4, 125 GeV) stau co- annihilation, right continent (3-8) 

(#24/7!", mo) = (0.58, 68 GeV) bulk (slepton-mediated), right continent (3.9) 



Model E is a typical stau co-annihilation scenario within mSUGRA, with a small mass difference 
m? 1 — mff = 4.7 GeV. It is beyond the reach of present LHC data, but will be eventually discovered 
in the jets+E 1 ™ 188 searches. The spin-independent nucleon-LSP cross-section is about a factor of 50 
smaller than what would be necessary to probe this model with the published XENON100 searches. 
Model F is also a stau co-annihilation model, but with a much more compressed spectrum. In this 
model rrifi — is less than 1 GeV. This small mass difference is needed in order for the co- 
annihilations to be efficient, given the larger LSP mass, and it means that the lighter stau will 
have only four-body decays f± — > u T V^iN\ and will be stable on time scales relevant for colliders. 
However, direct production of stau pairs is very small, and staus occur only rarely in gluino and 
squark decays in this model, so the initial discover at LHC will again come from jets+S™ lss searches. 
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FIG. 4: Comparison of superpartner and Higgs mass spectra for four selected model points E,F,G,H (top 
left, top right, bottom left, and bottom right, respectively) in the stau co-annihilation and bulk (slepton- 
mediated) dark matter regions in Figure [TJ 



The spin-independent nucleon-LSP cross-section is only a factor of 8 below the present XENON100 
limit for model F. Both models E and F have well above the LEP limit. 

Model G, in contrast, has asi about three orders of magnitude below the XENON100 limit, so 
dark matter direct detection will probably be impossible in the foreseeable future. Otherwise, it 
is similar to the CMSSM model E, with the main qualitative differences being a Higgs mass that 
is closer to the LEP limit of 114 GeV and a somewhat larger mass difference — = 7.0 
GeV. Model H is a bulk region model, with slepton-mediated annihilations N\N\ — > e + e~ , 
and t + t~ mostly responsible for dark matter annihilation in the early universe. The lightest stau 
is at mfj = 102 GeV, and about 13 GeV heavier than the LSP. Despite the low LSP mass, this 
model has asi more than 2 orders of magnitude below the present XENON100 bounds. The LHC 
signatures of models G and H should be very similar to model E, and a careful measurement of 
the slepton-LSP mass differences will be crucial for inferring their dark matter properties. 



13 




FIG. 5: Maps of the /i parameter (upper plot) and the predicted thermal relic abundance of dark matter 
f^DM^ 2 (lower plot) in the #24, mo parameter space with fixed M3 = — Aq = 600 GeV, tan j3 = 10, and fi < 0. 
The correspondence of shaded regions to in the upper plot and to fluM h 2 in the lower plot are the same 
as in Figure [T] The small corner near 6*24/^ = —0.07 and mo = 150 GeV is disfavored by the anomalous 
magnetic moment of the muon. 

It is also interesting to consider the same slicing through parameter space but with ji < 0. The 
map of \fi\ and Q^Mh 2 for this case is shown in Figure The qualitative features for this case are 
very similar to the positive \i case. In a small corner of the allowed central continent for mo < 200 
GeV and #24/^ near —0.07, the anomalous magnetic moment of the muon constraint eliminates 
some models; the other indirect constraints from eqs. (|2.4p - (|2.8p do not play a role for these models. 
Note also that in Figure the h° resonance region now is more pronounced on the right continent 
(rather than the left as for fj, > 0), providing solutions with Odm h 2 < 0.13 near 624 /it = 0.64 that 
are continuously connected to the bulk region. 

Figure [T] shows the ratio of the spin-independent LSP-nucleon cross section to the current 
XENON100 limit for that LSP mass, for model points in Figure Q] with J7dm^ 2 = 0.11. Here the 
most important feature that is different from the [i > case is that the small- models that are 
least susceptible to direct detection at XENON100 for /i < are those on the main, CMSSM-like 
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FIG. 6: The ratio of the spin-independent LSP-nucleon cross section to the current XENON100 limit for 
the model points in Figure [5] (with fixed M3 = —Aq = 600 GeV, tan j3 = 10, and /i < and varying 024, mo) 
that have QnMh 2 = 0.11. Different symbols are used for the model points according to which dark matter 
annihilation channels are most important in the early universe. 



continent. The focus point /small-// models on the right, large-#24 continent are nominally above 
the XENON100 limit, although we reiterate that this cannot be considered an exclusion. 

We now consider the impact of large tan /3 on the allowed parameter space. To illustrate this, 
we consider in Figure [7] the allowed regions in the (624, Mq) plane for tan/3 = 45, with fixed 
M3 = — Aq = 600 GeV and fi > 0. Here we note that the requirement of a neutralino LSP 
constrains the scalar masses mo to be larger than the minimal values obtained for moderate tan /3. 
In the central continent of Figure [71 we note a channel of points with ^dm^ 2 < 0.09; this is the 

-resonance funnel, where dark matter annihilation is driven chiefly by s-channel pseudoscalar 
exchange, with ^dm^ 2 = 0.11 ±0.02 regions on either side. This produces an island, centered near 
f?24/ 7r = 0.05 and mo = 700 GeV and separate from the main CMSSM-connected continent, which 
is a confluence between the stau co-annihilation, low-//, and A funnel mechanisms for efficient 
dark matter annihilation. However, this island is disfavored at 95% level by the constraint from 
BR(i? s — > < 1.1 x 10~ 8 . More generally, this constraint eliminates all of the stau co- 

annihilation and bulk region models on the central and left continents. On the right continent, 
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FIG. 7: Maps of the fi parameter (upper plot) and the predicted thermal relic abundance of dark matter 
^DM^ 2 (lower plot) in the #24, mo parameter space with fixed M3 = — Aq = 600 GeV, tan/3 = 45, and 
/i > 0. The correspondence of shaded regions to [i in the upper plot and to Hdm^ 2 in the lower plot are the 
same as in Figures [T] and [5j The regions under the contour in the lower plot are disfavored by the limit on 



however, a large number of stau co-annihilation models on the small-mo edge do survive, with 
values of /j, ranging up to 700 GeV. The top edges of the central and right continents are small-/x 
regions, as before, while on the left edge of the right continent, near #24/^ = 0.36, both s-channel 
A exchange and a significant higgsino content of the LSP play a role in dark matter annihilation. 

The ratio a si / oxenonioo for these models with tan (3 = 45 are shown in Figure [HJ Again we 
note that models with small fi > are most challenged by the XENON100 limits on the central, 
CMSSM-like continent. The models on the left side of the right continent, where the pseudoscalar 
exchange assists the higgsino content of the LSP in dark matter annihilations, tend to be far 
below the nominal XENON100 bounds even with the standard interpretation of those limits. As 
before, models with wino co-annihilation, ^4° resonance, h° resonance, and Z resonance as the 
dark matter annihilation mechanisms tend to be unchallenged by present direct detection limits as 
well. Although we did not encounter models in which the correct relic abundance of dark matter 
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FIG. 8: The ratio of the spin-independent LSP-nucleon cross section to the current XENON100 limit for 
the model points in Figure [7] (with fixed M3 = — Aq = 600 GeV, tan (3 — 45, and fi > and varying 6*24, mo) 
that have fioM^ 2 = 0.11. Different symbols are used for the model points according to which dark matter 
annihilation channels are most important in the early universe. 



is brought about by stop or sbottom co-annihilations 42|, |43j, there is no general reason why they 
should not exist, in particular if the assumption of a common scalar squared mass uiq were relaxed. 
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IV. EXPLORATIONS WITH FIXED M x 



In this section, we consider alternative slices in parameter space with M\ held fixed. This 
implies that in most of the parameter space, the LSP has a mass that does not vary greatly within 
the graphs to be presented below. With fixed Ml, the parameterization of eqs. (|2.ip - (l2,3p can be 
written as 

Therefore, fixing Mi implies that M2 and M3 become very large when 624/^ approaches —1/4 
and 3/4, so we again choose those asymptote values as the boundaries of the range, and again the 
parameter space splits into three main continents. At 624/11 ~ —0.102, one has M2 approaching 
0, so models near this vertical line are excluded, providing an ocean between the two continents. 
Similarly, at 624/^ ~ 0.148, one has M3 approaching 0, so models near this vertical line are likewise 
excluded, providing another gap between the two continents of viable parameter space on either 
side. The boundaries of the continents on either side of the M3 = line are actually determined 
by the fact that small IM3/M2I leads to small \n\ 2 , as can be seen from eqs. (jl.ip and (|1.2j) . 

We consider first models with fixed Mi = — Aq = 500 GeV, tan (3 = 10, and \x > 0. The maps of 
/i and Odm^- 2 for these models are shown in Figure EJ The allowed regions are qualitatively similar 
to the case discussed above for Figure [U except that in this case there is a region in which light top 
squarks mediate the dark matter annihilation, for 0.055 < 624/ V < 0.08 and 100 GeV < tjiq < 600 



GeV, as in ref. 30]. This region is now likely to be ruled out^ by LHC searches because the gluino 
is lighter than 660 GeV here, and the average squark mass is not much higher, although there is no 
experimental search specifically sensitive to this particular version of compressed super symmetry. 
More generally, the regions of the plane that the current LHC published searches are sensitive to 
are roughly indicated by the mq = 800, 1000 GeV contours shown in the lower plot of Figure EE 
Other regions have heavier squarks and/or gluinos. Note in particular that the left side of the 
right continent has gluino masses above 1.1 TeV, and increasing as one moves to the right. The 
region from 0.11 < #24/71" < 0.33 where one might have expected to find lighter gluinos was already 
excluded because —m/jj is negative at the TeV scale, precluding electroweak symmetry breaking. 
Because the gluino mass increases as one moves away from the center of the plots, jj, also increases, 
and unlike the fixed- M3 plots of the previous section, there are now models with [i well above 1 
TeV. This illustrates the fine-tuning price of moving to gluino and squark masses far above the 
LHC-accessible regions. 

Also indicated in Figure [9] are regions that are disfavored by the BR(6 — » 57) and limits. 
Recall, however, that the first of these can be evaded by intrinsically supersymmetry-breaking 
flavor violation. Moreover, it is now weaker than the probable direct reach of LHC in this particular 
parameter space slice. Note that these constraints eliminate the stau co-annihilation region with 
smaller fi on the CMSSM continent, but leave regions with fi > 450 GeV on the right continent. The 



However, similar models with non-universal scalar masses, in particular with larger gluino and first- and second- 
family squark masses, could evade these LHC searches. The LHC signatures for this scenario would include, besides 
the usual jets plus E ^p lss events, events with same-sign top quarks [3] and eventually striking resonant diphoton 
decays of stoponium [45l. |46| . a bound state of top squarks. Stoponium is long-lived enough to form bound states 
in these models because of the absence of two-body flavor-preserving decays of the top squark. 
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FIG. 9: Maps of the [i parameter (upper plot) and the predicted thermal relic abundance of dark matter 
^dmA 2 (lower plot) in the #24, mo parameter space with fixed Mi = — Aq = 500 GeV, tan/3 = 10, and 
/i > 0. In the upper plot, the dark (black) regions closest to the center represent // < 300 GeV, and 
successive regions away from the center correspond to 300 GeV < [i < 400 GeV (brown) and 400 GeV 
< fi < 500 GeV (red) and so on up to the last region (cyan) which corresponds to fi > 1000 GeV. Regions 
left blank do not have viable electroweak symmetry breaking, do not have a neutralino as the LSP, or have 
superpartners that are too light, as described in the text. In the lower plot, the thin dark regions (black) 
correspond to the observed range 0.09 < f^DM^ 2 < 0.13. The large interior regions (yellow) correspond 
to fiDM^ 2 > 0.13, while the darker shaded exterior region (green) has f^DM^ 2 < 0.09. Also indicated are 
contours for m h o < 113 GeV, for BR(fo — > s-f) from eq. (|2.7I) . and for to^ = 800, 1000 GeV corresponding 
roughly to the present LHC data reach, with regions closer to the center disfavored in each case. 

focus-point /small- fi regions on the top edges of the continents are untouched by these constraints, 
and we note that, as before, moving from the CMSSM focus point case to #24 > reduces the 
values of ran required for small /U and favorable dark matter, thus arguably substantially decreasing 
the fine-tuning cost. However, in the fixed M\ = 500 GeV parameter space, the constraint still 
imposes the requirement that mo is well above 1 TeV. Taking larger M\ would ameliorate this, as 
should be clear from the fixed M3 = 600 GeV cases studied above. 
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FIG. 10: The ratio ol the spin-independent LSP-nucleon cross section to the current XENON100 limit for 
the model points in Figure |9] (with fixed Mi = —Aq = 500 GeV, tan (3 = 10, and fi > and varying 6*24, mo) 
that have r2oM^ 2 = 0.11. Different symbols are used for the model points according to which dark matter 
annihilation channels are most important in the early universe. 



Figure [10] shows the ratio of the spin-independent LSP-nucleon cross section to the current 
XENON100 limit for the model points in Figure [9l For fixed Mi, we see the well-known fact that 
the correlation between the higgsino content of the LSP, inversely proportional to \f/,/Mi\ for these 
models, and the spin-independent cross-section is quite robust, with models having \i < 110 GeV 
strongly challenged within the standard interpretation of the LSP as dark matter (but not ruled 
out, as noted above). This includes especially the small-// models with smaller mo (i.e., the models 
near 624/^ = 0.07 and 0.37. Conversely, stau co-annihilation and wino co-annihilation models not 
otherwise ruled out are fine, and models with // > 300 are two orders of magnitude away from 
the sensitivity in direct dark matter experiments needed to challenge them, within this slice of 
parameter space. 

As a final exploration, we consider models with fixed M\ = —Aq = 900 GeV and tan/3 = 45, 
with [i > 0. The maps of [i and S^dm^ 2 are shown in Figure [TT1 as before. Notable features include 
the facts that regions with small /i are much thinner within this parameter space, and there are 
no viable regions with mo less than about 500 GeV. As in the case of fixed M3 with large tan/3 
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FIG. 11: Maps of the /i parameter (upper plot) and the predicted thermal relic abundance of dark matter 
Udm'i 2 (lower plot) in the #24, m parameter space with fixed Mi = —A = 900 GeV, tan/? = 45, and 
/i > 0. The correspondence of shaded regions to /i values in the upper plot and to f^DM^ 2 in the lower plot 
are the same as in Figure [9] Also indicated are contours as described in the text for the branching ratios of 
B s — > and B — > tv and b — > sj. 



in the previous section, there is a confluence island, on the shores of which ^4° resonance, stau 
co-annihilation, and small-// play a role in reducing the dark matter density. However, again the 
B s — > fi + n~ constraint strongly disfavors this island, and the smaller fi part of it is also disfavored 
by B — > tv and b — > 57. On the right continent, there is another large island nearly split off from 
the main continent by a channel consisting of an ^4° resonance funnel region. On the top edge of 
this island, pseudoscalar exchange plays the main role in the annihilation of dark matter, while on 
the lower edge it is mainly stau co-annihilation. The entire Qdm^ 2 < 0.13 part of this continent 
is free from the indirect constraints eq. (|2.4p - (|2.8p . It is also not challenged by dark matter direct 
detection, except perhaps for the top edge with small /x, as can be seen in Figure [T2l This is more 
generally true of all models within this slice of parameter space with /i > 180 GeV. 
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FIG. 12: The ratio of the spin-independent LSP-nucleon cross section to the current XENON100 limit for 
the model points in Figure [TT] (with fixed Mi = — Aq = 900 GeV, tan/3 = 45, n > and varying 6*24, mo) 
that have f^DM^ 2 = 0.11. Different symbols are used for the model points according to which dark matter 
annihilation channels are most important in the early universe. 



V. OUTLOOK 



As LHC searches continue to push the bounds on superpartner masses higher, there appears 
to be growing tension between the weak scale mz and the dimensionful supersymmetry-breaking 
terms that determine it from the effective potential. Non-universal gaugino masses provide a 
way of ameliorating this fine-tuning issue. In this paper we showed that there are large regions of 
parameter space in which as determined by the supersymmetry-breaking Lagrangian is naturally 
small, while scalar masses can be much less than the several TeV scale found in the traditional 
focus-point region in the MSSM. These models do evade the present bounds set by the LHC, but 
are now being confronted by dark matter direct detection experiments. However, exclusion is not 
possible with present XENON100 bounds for most of the cases we looked at, even if .R-parity is 
conserved and the lightest neutralino is really the dark matter. As is well-known, since the spin- 
independent cross section increases with higgsino content of the LSP, which in turn increases as 
decreases, future dark matter direct detection searches will probe the more natural models first. 
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In this paper, we also showed how non-universal gaugino masses can provide for a bulk region 
in which f-channel slepton exchange (without significant co-annihilations) plays the dominant role 
in annihilating the dark matter down to acceptable levels in the early universe. This is interesting 
because this region is severely endangered or perhaps extinct within the CMSSM after LHC and 
LEP bounds are considered. We further mapped out the regions in which stau co-annihilations and 
Higgs resonance play the most important role in dark matter annihilations; these regions differ in 
interesting ways from the corresponding CMSSM scenarios, and are not always continuously con- 
nected to them. In order to keep our st udy finite, we did not consider scalar mass non-universality 



here, but other works (see for example [22J and references therein) have showed that non- universal 
Higgs mass models have large variations in the dark matter relic density and spin-independent 
cross-sections. If supersymmetry is discovered at the LHC and in direct dark matter detection 
experiments, it will be a challenge to understand to what extent non-universalities in the soft su- 
persymmetry breaking parameters play a role in determining the phenomenological collider and 
dark matter properties of the theory. 

As this paper was being prepared, the ATLAS and CMS detector collaborations announced hints 



47j of a Higgs scalar boson signal in the vicinity of Mh = 125 GeV. If these hints are confirmed, 
the impact on the models considered here would be profound, as most of the parameter space 
cannot accommodate such large Mh values. If one assigns a combined 3 GeV uncertainty (note 
that this is a somewhat arbitrary estimate, including both the theoretical uncertainty of about 2 
GeV [48] and the quite preliminary LHC uncertainty) to Mh, and therefore considers regions with 
predicted Mh > 122 GeV to be acceptable, then among the models we considered only those with 
moderately negative 024 (on the main continent) and large too are viable. From the point of view 
of the putative Mh = 125 GeV signal, the variation on the usual CMSSM focus point region with 
negative 024 is therefore preferred. There are several other ways to increase the prediction for Mh, 
including by raising all of the soft masses (by raising M3 in our framework) and by increasing 
tan /?. However, among the models considered here, the heavy scalar solutions with negative 024 
would seem to be the easiest to make consistent with Mh ~ 125 GeV. Other mechanisms, such as 



adding new vector-like supermultiplets 49|]- 53J, would affect our results in qualitative ways. 
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